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Abstract 

We perform a non-perturbative chiral study of the masses of the Hghtest pseudoscalar mesons. In the 
calculation of the self-energies we employ the S-wave meson-meson amplitudes taken from Unitary 
Chiral Perturbation Theory (UCHPT) that include the lightest nonet of scalar resonances. Values for 
the bare masses of pions and kaons are obtained, as well as an estimate of the mass of the rjs- The 
former are found to dominate the physical pseudoscalar masses. We then match to the self-energies 
from Chiral Perturbation Theory (CHPT) to ©(p*), and a robust relation between several 0{p*) 
CHPT counterterms is obtained. We also resum higher orders from our calculated self-energies. By 
taking into account values determined from previous chiral phenomenological studies of m^/m and 
3Li + Ll, we determine a tighter region of favoured values for the 0{p'^) CHPT counterterms 2Lg — L\ 
and 2Lg — Lg. This determination perfectly overlaps with the recent determinations to 0{p^) in CHPT. 
We warn about a likely reduction in the value of nis/m by higher loop diagrams and that this is not 
systematically accounted for by present lattice extrapolations. We also provide a favoured interval of 
values for rris/m and 3Lj + Lg. 
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1 Introduction 



Chiral Perturbation Theory (CHPT) is the effective quantum field theory of strong interactions at low 
energies [ll[2l[3]. Three flavour CHPT, including strangeness, has already reached a very sophisticated 
stage with present calculations at the level of next-to-next-to-leading order (NNLO) or 0{p^) [4l[5l[6]. 
In ref.[71[8] the masses and decay constants of the lightest octet of pseudoscalar mesons are worked out 
at to 0{p^). In this reference it is shown that the 0{p^) two loop contributions to the self-energies of the 
pseudoscalars are much larger than the full 0{p'^). In addition, the 0{p^) tree level contributions to the 
kaon and eta masses, with the 0{p^) chiral counterterms estimated by resonance saturation, are much 
larger in modulus than the 0{p^) pure loop ones and with opposite sign. As stated in ref.[8] the 0{p^) 
counterterms "seem severely overestimated" due to the complicated nature of the scalar sector. Hence, 
this reference rises the interesting question of how one can improve the calculation of the contributions 
from the lightest scalar resonances. It is well known that the scalar sector is a source of large higher order 
corrections to the CHPT series. As a relevant example, let us quote the / = tttt scattering length, Og, 
which receives large higher order corrections (a 40% correction with respect to the lowest order value). 
This should be compared with the non-resonant 1 = 2 S-wave scattering length, Og, which has negligible 
higher order corrections [9l [T0\ [TT] . Even larger higher order chiral corrections due to the 1 = S-wave 
final state interactions happen for the lowest order prediction of T(r] Svr), more than a factor of 2 
[12^ I13j . Another good example is the process 77 tt^tt^ for which a two-loop CHPT calculation is 
needed [14:\ I15j to improve the comparison with dispersive theory calculations |16l 1171 [T8]. connection 
with these enhancements one has the generation of the lightest scalar resonances a, /o(980), ao(980) and 
K [IHl Uni [2TI [22] because of the self-interactions among the pseudoscalars in S-wave. 

The Feynman-Hellman theorem |23] relates dM'^ /dnig, where M is a hadronic mass and niq the q 
quark mass, with the qq scalar form factor of this hadron at zero four-momentum transfer, the so called 
sigma-terms [24]. Of special relevance are those in connection with the lightest hadronic states, tt, K, r] for 
mesons and the nucleon for baryons. The large impact on the evaluation of the pion-nucleon sigma-term 
of the I = J = TTTT amplitude, due to the ^-channel tttt exchange, is remarked in ref.[25]. On the other 
hand, it was stressed in refs. j26l [271 I^S j I29j that the scalar form factor of the pion receives large unitarity 
corrections so that the one-loop CHPT approximation becomes inaccurate at a surprisingly low energy. 
Here, the so called infrared singularities are enhanced in comparison e.g. with the vector form factor 
|27[l20j. In ref.[29] the scalar KK form factors were first evaluated at the one loop level in CHPT and also 
employed in UCHPT. The corrections to these form factors because of the final state interactions of the 
/ = J = TTvr and KK coupled channels were huge, see figs. 5 and 6 of ref.[29j. Since this happens to the 
scalar form factors one should expect a similar large impact of the I = J = pseudoscalar-pseudoscalar 
amplitudes in evaluating the pseudoscalar masses, as both quantities are related by the Feynman-Hellman 
theorem. We offer here an estimation of such effects, that includes the exchanges of the lightest scalar 
resonances a, k, /o(980) and ao(980). 

Recently, precise calculations of the pseudoscalar masses in lattice QCD with three dynamical fermions 
are available [30[ \3T\ [32] . In these evaluations lattice results are taken down to the values of the lightest 
quark masses, and m^, and extrapolated to the continuum by employing perturbative SU(3) Staggered 
CHPT calculations [33l 131] • Given the large size of the pure loop contributions at 0{p^) estimated in 
ref.[8], and since they were neglected in refs. [311 [32], one should consider the possibility whether these large 
contributions are buried in the present values for quark masses and low energy parameters determined 
by this Collaboration. 

The manuscript is organized as follows. Section [2] is dedicated to the development of the formalism for 
the calculation of the self-energies. In section [3] the results without matching with CHPT at O(p^) are 
given, this is what we call the full dynamical self-energies. The 0{p^) CHPT self-energies are introduced 
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and supplemented with higher order corrections from our formahsm in section HI We also derive in this 
section a robust relation between -^^(4,6) = 2Lg — L4 and ^(5,3) = ^Lg — -L5 and determine our interval 
of values for these counterterms. It turns out that this region remarkably overlaps with the values 
determined in refs.[35', '7] from 0{p^) fits to data. We discuss in section [5] about likely effects from higher 
loop diagrams and determine the values of ms/m, -^^(7,8) = 3L7 + Lg, and self-energies from our favoured 
region of low energy chiral counterterms. Some conclusions are collected in section [H 



2 Formalism 

Our starting point is the S-wave meson-meson partial waves both for the resonant / = 0, 1 and 1/2, 
as well as for the much smaller and non-resonant ones, 7 = 3/2 and 2. For the former set we take the 
amplitudes of ref. [20] . The interaction kernel employed in this reference comprises the lowest order CHPT 
amplitudes together with the s-channel exchange of scalar resonances in a chiral symmetric invariant way 
from ref.^36j. These tree level resonances constitute an octet with mass around 1.4 GeV and a singlet 
around 1 GeV. On the other hand, the coupled channels are vrvr, KK and rjrj for / = 0, vrr/ and KK for 
1 = 1 and Kir and Krj for / = 1/2. In addition to the resonances explicitly included at tree level, related 
to the physical ones around 1.4 GeV like e.g. the Kq{14:10) or the ao(1450), the approach also generates 
dynamically the resonances a or /o(600), k, ao(980) and the main contribution to the /o(980). These 
resonances are generated even when no explicit resonances at tree level are included. The basic point 
in UCHPT is to resum the right hand or unitarity cut to all orders, the source of the large corrections 
produced by the S-wave meson-meson interactions, and perform a chiral expansion of the rest, the so 
called interaction kernel [371 138^ [20j. For the non-resonant isospins 3/2 and 2, not given in ref. [20], the 
kernels correspond to the O(p^) CHPT amplitudes. In fig{T]we show the reproduction of scattering data 
that is achieved by our S-wave amplitudes. 

We only consider the S-waves since the lightest pseudoscalar self-energies are expected to be domi- 
nated by low energy physics where the P-waves are kinematically suppressed by small factors of three- 
momentum squared. This is true even in the calculations of loops since, after proper renormalization, 
the typical momentum in virtual processes will be bounded and of the same order as the on-shell values. 
Indeed, by a comparison of the order of magnitude of our results with those obtained by regularizing the 
loops with a three-momentum cut-off, the latter must be around 0.3 — 0.4 GeV. At such energies only 
S-waves matters in very good approximation, which is also clear experimentally. Another signal of the 
suppression of P-wave dynamics in the self-energies is the fact that it only starts to contribute at the 
two chiral loop level. In addition, we recall the discussion in section [T] concerning the large corrections 
to the chiral series typically induced by the resonant S-wave channels and the large 0{p^) counterterms 
calculated from resonance exchanges in ref. [8]. 

There is recently both experimental and theoretical overwhelming evidence of the existence of the 
I = a [39l[l9l[^|40l[4ll|42] and/ = 1/2 k [39l[20l[40l[2ll|42l|43] resonances. One also has the 
1 = resonance /o(980) and the 1 = 1 ao(980). These resonances have been considered quite often in 
the literature to form the lightest nonet of scalar resonances [441 l20l HOl 145] . though there is not still 
consensus on this issue. Since these resonances are the lightest scalar ones it is an interesting question to 
calculate their contributions to the pseudoscalar masses. In connection with this, we plot in figl2| those 
contributions to the pseudoscalar self-energies^^at follow from the exchange of scalar resonances R at 
the level of one loop with dressed propagators 
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*'^There is a similar contribution to that of figOJa with the resonance leg coupling directly to the vacuum without the 
tadpole in the upper extreme for the preexisting octet of scalar resonances at 1.4 GeV. These contributions were included 
e.g. in ref. [21 1 [46] and their contributions to the chiral series for evaluating self-energies are negligible as they only contribute 
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Figure 1: Different phase-sliifts, inelasticities and mass distributions obtained with tlic unitarized amplitudes. 
The experimental data are from refs. [171llllli[5ni[ni[51[51[Ml[S3[51[571[Ml[5S] 



Within UCHPT [5U] the resonances /o(980), ao(980), a, k and an octet around 1.4 GeV appear as 
poles of the corresponding S-wave amplitudes once the right hand cut is resummed according to figl3^. 
In this way figl2^ is included in the evaluation of figl3}3 and, analogously, figl5)o is a contribution in 
figSb- As remarked above, the lightest scalar resonances are generated through the self-interactions of 
the pseudoscalars \19\ [20] , with the kernel given by Chiral Symmetry [H [21 [3] . These scalar resonances in 
figs[2] better correspond to two pseudoscalars strongly correlated by the strong self- interactions. Hence, 
in order to avoid double-counting one cannot talk independently of the contributions of the unitarity 
loops and the light scalar resonances. Both must be calculated at once. This is why a scheme like that in 
figs [3] and [His the appropriate one. For the resonances higher in mass, namely, the octet around 1.4 GeV, 
the situation is quite different as they are preexisting resonances that are dressed by the interaction with 
the pseudoscalars. In this way, tree level exchanges and pseudoscalar loops can be differentiated in a 
perturbative way of thinking. 

For the calculation of figlSb we proceed by analogy with figi2^. Note that in ref.[6p| it was shown 

to 2Lg — LI and 21/5 — L4, modulo small higher order contributions, with a vanishing result. 
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Figure 2: Fig. 2a represents the s— channel exchange of scalar resonances R in terms of the intermediate Q and 
initial P pseudoscalars. A sum over all possible intermediate Q and R should be understood. The fig. 2b corresponds 
to the crossed exchange of the / = scalar resonances R. 



that the pion and Kt: scalar form factors are dominated by the Laurent poles of the a and k resonances, 
respectively, for energies below the next resonance. In the case of the / = and 1/2 S-waves it was 
similarly shown that they are given by the previous poles plus a background, which mainly consists of 
just a constant which would give rise to a term analogous to figlDa or ^p. It follows then that figl2^ is 
a close analogy to figl3)3. After performing the integration over the null component of q in the loop, see 
figl2^, one has two type of cuts. The first corresponds to having the pseudoscalar Q on-shell, while for 
the second the state on-shell is -R, running with opposite sense to that of R. The former involves the 
PQ scalar interaction mediated by the scalar resonance R. This is the contribution we are seeking. The 
latter corresponds to the PR pseudoscalar interaction mediated by the pseudoscalar Q which is wildly 
off-shell and suppressed because the resonance energy is —En. 
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Figure 3: Fig. 3a represents the S-wave amplitude PQ PQ and fig. 3b is the diagram for the calculation of the 
self-energy of the pseudoscalar P due to the intermediate pseudoscalar Q. 
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Figure 4: Tadpole diagram for the self-energy of the pseudoscalar P. 
We then isolate from figs 12^ and[3j) the following scalar contribution to the self-energies of the pseu- 
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doscalars, 



yU 



E 



(2.1) 



with Q = {tt+ ,iT- ,iT° , K+ , K~ , , K ,7]} and si = {Mp - EQ(k)f - k"^ for p = (Mp,0). The previous 
equation is invariant under Lorentz transformation as must be the case for a self-energy because a meson- 
meson scattering amphtude is a Lorentz scalar. The sum in eq. (12.ip is over all the species of pseudoscalars 
because all of them give rise to different PQ^ — > PQi amplitudes. 

As stated above, one also has the corresponding diagrams figsi2)3 andHja. Notice that in figl3]we were 
driven to consider the amplitudes PQ^ — > PQi, but one has to take into account the t— crossed process 
PP — > QiQi as well. We can evaluate the diagram in figlD3 by performing the contour integration for q^. 
Here one only has the singularity due to the Q pseudoscalar and, 



-<tad 



A, 



Q 



Q=ttO,tt+,K+,KO,' 



(27r)32i?Q(k) PP-QQ 



(2.2) 



In Sp*^ one sums only over particles (without including the sum over antiparticles) since in the final 
state both particles and antiparticles happen and otherwise one is double-counting. Finally, Aq = 1 for 
7r+, and and Aq = 1/2 for vr*^ and rj, because they are their own antiparticles. These warnings 
also happen in the standard calculations of tadpoles with the use of Feynman rules because of the same 
reasons. The self-energy of the pseudoscalar P is then given by the sum, 



lup — lup + lup , 

In terms of the amplitudes for the different channels we have for the vr, K and r] self-energies, 



(2.3) 
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:iT^-i^^{s) + AT;- 



1=1/2 



^1=0 



V3T^-\^{0) - 2r4-_^^^(0) + T^-^rj^(O) 



(2.4) 



Notice that we are employing just the S-wave contribution to the full strong amplitude Tjj, and this is 
why we have kept only the s-Mandelstam variable in the argument of Tij{si), dropping the t dependence 
in eq. (j2.4p . This is what automatically occurs when considering the intermediate scalar resonances in 
figsi2^, b. The Clebsch-Gordan coefficients in front of the S-wave amplitudes in the previous equations 
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are obtained from a standard isospin analysis. As an example, for the tt"^ self-energy one has the S-wave 
amplitudes, 

-'■TT+TT" — J-TT+TT" 7r+ 7r+ — >7r+ 7r+ TT+TT*^— J-TT+TT*^ 1 n TTTT — i-TTTr"'" O TTTT — i-TTTT I TTTT — >7r7r' o TT 



g ^ TTTT — i-TTTT ' ^ TTTT — i-TTTT / ' TTTT — ^TTTT ' ^ '^'^ — >7T7T 

To finally obtain the coefficients in Syr one has to multiply by 2 the T^^_^^^ S-waves in the last expression 
because in refs.[l9l [20] one uses the / = 0, 2, tttt states (and also the I = rjrj one) normalized to 1/2. 
In this way, these symmetric states under the exchange of the two pions can be treated as if the latter 
were distinguishable. For more details see refs.|19t I20j. 

The integral in eq. ()2.ip for Sp is ultraviolet divergent since for large three-momentum the measure 
grows as and TpQ^pQ{si) tends to 1/logsi or constant for |si| oo. This, together with the factor 
Eq{]<.) in the denominator, gives rise to a quadratic divergence. For Sp'^, eq. ()2.2p . one has the same type 
of divergence as Tp-p^QQ is evaluated at si = and it is a constant. 

Q 



p P 

Figure 5: The one loop contribution of figslSb andHb at 0{p'^). This is the only loop appearing at C(p^) CHPT. 

A general coupled channel S-wave T-matrix has the structure |20j . 

T = [/ + /C-G]"^-/C , (2.6) 

with /C the interaction kernel and G corresponds to the unitarity bubbles shown in figl31 For more details 
see ref . |20j . The previous equation has the expansion, 

r=[/-/C-G + /C-G-/C + ...]-/C. (2.7) 

We notice that when T = IC2 (the lowest order CHPT amplitudes), one has no unitarity bubbles in 
figEb and figSb, and then both figures give the same diagram. This only happens for IC2 at lowest 
order in CHPT. The local term IC2 gives rise to the one loop diagram shown in figlHl the only one that 
appears at next-to-leading order (NLO) in CHPT (at this order the rest of contributions come from local 
counterterms [11 [3]). In order to avoid this double-counting we then subtract the contribution in figl5]to 
the integral in eq. (j2.ip . fig|3)3. This is done explicitly below in eq. (j2.27p . 

2.1 Regular izat ion and renormalization 

We now proceed with the evaluation of the integral in eq. (121^ for S^. The integral in eq. ()2.2p for Sp 
is a particular case of the previous one. Firstly, we perform the angular integration and change to the 
energy as the integration variable. 
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We now introduce the new variable 

t = Mq/Eq . (2.9) 

As the resulting integral diverges quadratically for t — > 0, we change the lower limit to e, positive and 
small, taking at the end the limit e — > 0"*". 

^ ' dt^L^Tp^ {s,) , (2.10) 



with si = M|, + - IMpMq/t. Notice that since e = Mq/A, with A ^ Mq being the upper limit in 
the three- momentum of eq. (12.8p . the previous limit is equivalent to A ^ +oo. In performing this limit 
we must identify those terms that scales as 1 /e"^ , 1/e (equivalently as A^ , A) and remove them from the 
result. In order to remove infinities one must work out algebraic expressions for the previous integral that 
explicitly show these diverging powers of e. The problem of accomplishing this aim is the complicated 
matrix expression for the strong amplitude that prevents us from giving a close expression for the integral 
in eq. (l2.10p . However, we are interested in isolating algebraically the dependence in e, and this can be 
done by performing a dispersion relation representation for TpQ^pQ^si) in the physical Riemann sheet. 

The Tpq_^pq{si) amplitude from ref.[20] has only one cut, the right or unitarity one, and tends to 
1/logsi or constant when si — > oo. In order to perform the dispersion relation we take as integration 
contour a circle in the infinity deformed to engulf the real axis along the right hand cut, Sfh < si < oo, 
with sth the lightest meson-meson threshold with the PQ S-wave quantum numbers. We take one 
subtraction because of the aforementioned behaviour of Tjj(si) at infinity. We then have: 

Tij[si) = Tij{s2) + — - — / ds — — — — + > — 7j — -7j , (2.11) 



TT 



(S'-Sl)(s'-S2) fri (sf - Sl)isf - S2) 



with S2 the point where the subtraction has been performed. In addition, N is the number of poles of 
Tij{s) in the physical Riemann sheet and R^^^ is the residue of the £th pole at the position sf. Because 
of the Schwartz reflection principle any partial wave satisfies that Tij{s*) = Tjj(s)*, so that poles always 
appear as pairs in relative complex conjugate positions with complex conjugate residua. Taking then S2 
real in eq. (|2.1ip the contribution from the sum of poles is real for si along the real s-axis. 

A physical partial wave should not have any pole on the physical sheet with non- vanishing imaginary 
part. But notice that the lowest order CHPT S-wave amplitudes in / = 1/2 {Ktt Kit, Kir Kt] 
and Kt] — > Kt]) have each of them one pole at s = due to the kinematical reason of having different 
masses, see ref.[20] for explicit expressions of those amplitudes. This pole at s = for the interaction 
kernel drives the appearance of a pole in the final Ktt — > Kit S-wave for a real s value, typically at 
around s = —0.25 GeV^, not far away from its lowest order value. This poles gives rise to a non- 
negligible effect in the calculation of the vr self-energy. There are also other poles with non-zero imaginary 
part. Reassuringly, these poles are always far away of the physical axis and their contributions to the 
pseudoscalar self-energies are negligible, as we have numerically checked. 

Let us proceed with the evaluation of the renormalized value of the divergent integral eq. (j2.10p making 
use of the representation given in eq. (j2.1ip for the S-wave T-matrix elements. Both the contributions 
from the unitarity cut and poles involve the same sort of integral, 

J{s,S2)= I dt^^y-^. (2.12) 
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For the sum over the poles in eq. (|2.1ip the variable s' must be replaced by sf. For convenience, new 
variables a and b are introduced, 

a = u — ,b = u , (2.13) 

2MpMq ' 2MpMq ' ^ ^ 

with I' = {Ml + M^)/{2MpMq). Notice that one can also write l/t = - si/2MpMq. Then, 
J{s',S2) = - dt——^ — T- = - — + {a-b) dt- 



1-bt J, ' 'J, 

+ (a - b)b / dt— — - - (a - b)b / . (2.14) 

Je t Jet- 

The first integral is elementary, although its renormalized value requires some comments. 

dt^^ = -^-^ + 2 log ^ . (2.15) 



£ 



Now, when sending e ^ O"*" both terms on the r.h.s. diverge. Performing a power series around e = 
one has — l/2e^ + (1 — 21oge/2)/4 + Oie^^ The first term in this sum diverges like and should 
be reabsorbed by appropriate chiral counterterms. The terms with positive powers of e, not explicitly 
shown, vanish for e — > 0. Then, the remaining contribution is, 

7 1 - 21og^ =71 - 21og^ - o log^ , (2.16) 



4V 2Ay 4V 2/xy 2°A 

where the renormalization scale \x is introduced. In sending A +00 (e — > 0"^), again the term log^/A 
diverges and should be reabsorbed by the appropriate counterterms. Then one has: 

where the symbol =, used for expressing our renormalized integrals, means equal up to terms divergent 
as e~^, or loge, when e — > O"*". When giving our results we shall vary \jl between 0.5 and 1.2 GeV, 
providing in this way an estimate of possible constants of order one that may change depending on the 
renormalization scheme chosen. E.g. the difference between our procedure and the results in the modified 
MS scheme in dimensional regularization used in CHPT [21 [3]. This variation in [i will be a source of 
uncertainty in our results and it will be consider in the error analyses. In ref.[39j the unitarity bubble 
G in eq. ()2.6p was also calculated within a cut-off scheme. If one sends this cut-off to infinity only a 
logarithmic divergence remains, as we have here. By varying within the range mentioned above, one is 
accounting for finite constant terms, like the subtraction constant used in G. 

The second and third integrals on the right hand side of eq. (j2.14p can be treated similarly. For the 
last one, 

dt ^^~*^ = Vl - t2 - c arcsini + (1 - c^) [ , (2.18) 

t-c J (t- c)^l-t'^ 

with c = 1/6 and 



lib) = / , = ^) — , log , 6 > 1 

^' ' {t-i/b)VT^ Vi - i/fe' ^ + v& 
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... 1 ^ 1 + - 1/62 

zx) , log ^ ^ ,' — , 6 < -1 , 

^^^^^ _ ^ (I + arcsinft) , < 6 < 1 , 

iv) , ^ f 5 + arcsin h] , - 1 < 6 < , (2.19) 
' V^l/62 - 1 V2 y ' > V y 

Adding all the integrals that contribute in eq. (|2.14p . 

J{s', S2) = ^ + ^ log - + (a - 6)6 log - - (a - 6)6(1 - 1/62) J^^^ _ ^2.20) 

Now, in the expressions given for / with |6| > 1 in eq. (|2.19p . the first two lines, one has log(e + 1/|6|) = 
loge + log(l + l/e|6|). This last term gives rise to an infinite series of divergent terms in powers of 
1/(|6| e)" for |6| > 1, equivalently in powers of (A/6)"', which are removed in the regularization process. 
Substituting also loge by log Mq/ fj,, introducing the renormalization scale fi, as explained before, then 

J{S',S2) = 

i) J+{a, 6) = i + i log ^ + (a - 6)6 log 2 - (a - 6)6(1 - ^1 - 1/62) log 



42° Mq^ '° ^ ' ' " ' ' ° fi 
-{a - 6)6^1 - 1/62 iog(i + ^1 - 1/62) , |6| > i , (2.21) 
1 1 2/i 2/i 



ii) J -{a, 6) = - + - log ^ + (a - 6)61og ^ - (a - 6) Vl - 62 (| + arcsin 6) , |6| < 1 . 
The dispersive integral contribution to the self-energy of the pseudoscalar P by the intermediate Q is: 



i) If (Mp - Mq)2 > Sth 



Jstu {S' - S2) vr J(^Mp-Mq)-^ {s' - 82) 

1 /•+~ , ,Imr,,(s') , , , , 

TT J{^Mp+MqY (■S - S2) 
ii) If {Mp - Mg f < Sth 

i/ + i / dr^^^^,U{s\s2) (2.22) 



vr 



(s' - S2) vr J(Mp+MQy {s' - S2) 



Notice that {Mp + Mq)"^ > sth- Let us remark that J+(s',S2) vanishes like — 1/s' for s' ^ 00, as it is 
clear from eq. (j2.2ip . This is why all the last integrals in the previous equations converge for s' +00. 
The subtraction constant in eq. (12.1ip . Tij{s2), contributes as. 



CI = T^J{S2) I di^^V- = -^^.(^2)7 ( 1 - 21og ^ ) , (2.23) 



— = -r,te)j(i-2iog-; 

from eq. (l2.17[) . The contribution of the sum of poles in eq. (|2.1ip is: 



e=i '^'^^ •^2j 
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We have only J+(a, b) since for this case always |6| > 1 as the poles happen far away in the complex plane 
or on the negative s-axis. If S2 is taken such that a = then J(a, b) will vanish like 0{b~'^) for \h\ oo, 
as follows from eq. (|2.2ip . One can take advantage of this fact since the pole contributions will tend to 
vanish because sf and the resulting \b\ are large. Thus, we fix S2 = sa = Mp + Mq in the following 
so that a = 0, eq. (|2.13p . This choice also makes that all the contributions in our final equation for the 
self-energies, eq. (l2.29p below, have natural size. 

One still has to remove the 0{p^) contribution to -fj^^ in order to avoid double-counting with the 
tadpole contribution, as explained above after eq. (j2.7p . We then evaluate eq. (j2.ip with Tpg^pg given 
by its 0{p'^) expression. A general CHPT S-wave amplitude at 0{p^) has the form. 



T^f{s) = A + Bs + -. (2.25) 



C 



s 



The last term in the sum only appears for the 1=1/2 Kir and Kr] coupled channel amplitudes involving 
two particles with different masses. Then, 



dt — — [A + Bsi{t) + 



L ' ' Slit), 

= {A + saB) - 2MpMqB dt^^^ + £ dt^^^^^ , (2.26) 

with f = l/u (f < 1 when C ^ 0) and sa = ^p^Mq. The renormalized second integral in this equation 
vanishes as can easily worked out. We then have. 



Irem " {sa) ^^^^ + " ^3 t_f 



= -r,f (.^)i(l-21og^) + £j+(0,l//) . (2.27) 
The other contribution to Sp is S^'^, ea.dOD. Its calculation, eq. (j2.15p . is straightforward and gives. 



e 



Hence, for eq. (|2.3p one has: 

So - - V — ^ (l^^ + I^^ + I^^ - I^Q") - V —^I^Q' (2 29) 

— ^2t:)2 V^'^bs^ ^disp^ ^pole ^remj (^27r)'^ ^'"^ l^-^yj 

with Q = {Tr°,7r+,7r~,K^,K'^,K+,K~,r]} and Q' = {n^ ,tt+ , K+ , ,7]} . Notice that after the sub- 
traction of YIq -^rem MQ/{27r)'^ to the previous equation, the only contribution at 0{p'^) of our result 

comes from lfJ^M^/{27rf when Tij{0) in eq. (j2.28p is substituted by tI^j\o), the corresponding low- 
est order CHPT amplitude. It is important to remark that for this 0{p'^) case we have checked that 
'YliQ' -^tad -^Q'/(2^)^ reproduces the CHPT infrared logarithms on the quark masses in CHPT at 0{p'^) 
[3]. This is one among the many explicit calculations in the literature showing that by sending A — > -|-oo 
and removing the power divergences (which necessarily disappears in dimensional regularization) one 
recovers the dimensional regularization results with the differences reabsorbed in the chiral counterterms 
from the tree level contributions |6H [621 [63] . 
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2.2 Mass equations 

We often remove the 0{p^) contribution to Sp and denote by Sp the 0{p^) and higher order contri- 
butions. The latter is obtained with the replacement of Tij(O) by Tij{0) — T>j (0) in eq. ()2.28p for the 
tadpole contribution. The O(p^) CHPT expressions for the pseudoscalar masses are, ref.[3]: 

= U + /u^ - -/i^ + 2mK3 + 

r2 2 



M'k = i 1 + 3^'7 + ("i + '^s)K^ + Ki \ , (2.30) 



9 o2 



where, 



= -^{2Ll-Ll) , K, = {2m + ms)—^{2U,-Ll) , 
^5 = K-m)^^||(3L, + Lg),;.Q = ^|flog^. (2.31) 



The symbol Mp denotes the bare mass of the pseudoscalar P corresponding to the lowest order CHPT 
result, 

M^=2Bom, Mk= Bo{7h + ms) , M„= -So(m + 2ms) , (2.32) 



and Bq measures the strength of the quark condensate (0|g(7|0) = —FqBq + 0{'mg) in the SU(3) chiral 

. o 2 

hmit. We will denote the 0{p^) CHPT self-energies by Sp^ = Mj,- Mp, eq.ig^D]). Notice that from 

o 2 

eq. ()2.32p always satisfies the Gell-Mann-Okubo mass relation. 



o 2 1 o 2 o 2 



M^= -(4 Mi^ - M^) . (2.33) 



On the other hand, the -L[ coefficients are the renormalized values of the low energy CHPT counterterms 
at 0(j)^) from ref.[3] within the MS — 1 scheme. 
The mass equation one has to solve is, 

2 2 

M| =Mp +^p{Mq; M|) , (2.34) 

o 2 

with Ep given in eq. (|2.29p . The dependence on the bare masses Mq originates from the one of the 
interaction kernel /C |20j . 

If we impose that our results incorporate the 0{p'^) CHPT self-energies and resum higher order con- 
tributions, the following equations arise. 

Ml =Mp +T.f{M% LI) + (Mq; ) . (2.35) 

Our analysis is performed in the isospin limit, niu = md = rh and no electromagnetic contributions. 
The symbols Mp above refer to the masses of the lightest pseudoscalars in the isospin limit. To determine 
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these masses we employ the Dashen theorem [M] which states that at the lowest nontrivial order in 
the neutral particles ir^ and do not receive electromagnetic contributions, while these are the same 
for and . In this way, 

= Ml, ~ (135 MeV)2 , M|- = (M^+ + M|o - MI+ + Mln)/2 ~ (495 MeV)^ , (2.36) 

are taken as the pion and kaon masses in the isospin limit, respectively. Violations of the previous 
relation for the pion mass are expected to be tiny, 0{{mu — nidf') [3], from the pure QCD side, and 
C'(e'^m^/(87r^/^)) from electromagnetic interactions for the vr'^ mass. More substantial seems to be the 
violations of the Dashen theorem for the kaon masses, of order e^Mj. [651 [Ml EI] . So that, if we write 
M]^ = [M^+ + M^o - (1 + 5D){Ml+ - M^o)] /2, with 5D measuring the deviation with respect to the 
Dashen theorem, then as a conservative estimate, 5D is expected to be in the range — 2 |30j from 
refs.[551 [57] . But even so, the changes in are at most a 0.5%, which can be discarded within the 
definitely larger uncertainties of our calculations, to be shown in the next sections. 



3 Dynamical self-energies 

o 2 

We restrict in this section to eq. (|2.34p and solve for M^^k (equivalently, for Bqiti and BquIs) and for 
M^. In this section the r] is actually the r]s since we are not including any source of tj — r]' mixing. The 
distinction between the r] and r/g is accomplished in the next section. 

In perturbative studies one solves iteratively eq. ()2.34p . In this way, for the first round, we will use 

o 2 

physical masses in all the arguments of the self-energies T,p{Mq', Mq) with the resulting values 

M^= 108 ± 4 , Mk= 422 ± 18 , M^g = 660 ± 17 MeV. (3.37) 
The errors are calculated by performing a Monte-Carlo sampling of the free parameters of the T-matrix 



of ref.[20^ and varying the renormalization scale // introduced above in the range between 0.5 and 
1.2 GeV. This will be the standard procedure to evaluate errors in our work and should be understood 
in the following. Performing one more iteration in eq. ()2.34p . i.e. using the previous bare masses in the 
appropriate arguments of the self-energies, it results 

M^= 158 ± 7 , Mk= 511 ± 12 , M^g = 627 ± 15 MeV. (3.38) 

o o 

We observe significant differences between eqs. (13.37j) and (13.38P for Mn and Mk, which is a clear 
indication that higher orders corrections, of 0{p^) and superior, are relevant for the calculation of the 
lightest pseudoscalar self-energies in SU(3). This was also observed in refs.[71[8], where the masses were 
calculated at 0{p^) in CHPT. 

Because of the large variation between the values in eqs. ()3.37p and (j3.38p . eq. (j2.34p should be solved 
exactly. This is done by minimizing the function, 

, m!+s.(Mq;M|) \ ( mj,+i:k{m%mI) 

S{M.,Mk) = I j^, 1 + 



2 2 

M^ +T.^{Mq;MI)-MI 
Ml 



(3.39) 



*■* We have made a new fit, including some more recent data, and tlie resulting fit has central values for the free parameters 
that are compatible with those of ref . [20) . within the shown errors in this reference. 
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This method is more robust than the iterative one and it can provide solutions even when the later 
does not converge. On the other hand, since we are now varying the bare masses employed in the kernels 
for the T-matrices, we refit the scattering data in fig(T]for every value of the given bare masses. Then, the 
T-matrix changes as a function of the values employed in each step. The following results are obtained, 

M^= 126 ± 4 , Mk= 476 ± 10 , M^g = 635 ± 15 MeV. (3.40) 

One observes significant variations between the exact result and the first and second iterated solutions, 
eq. (|3.37p and (j3.38p . respectively. This point is usually overlooked in the literature [3l[8]. Regarding the 
sizes of the self-energies we obtain 

° 2 

— ^ = 0.11 ± 0.06 , 0.08 ± 0.04 , 0.26 ± 0.06 , (3.41) 

for pions, kaons and etas, respectively. We then observe that the physical masses of the pseudoscalars 
are dominated by the bare masses and that the dynamical contributions because of their self-interactions 
(which include the exchanges of the scalar resonances cr, k, /o(980) and ao(980), as well as of the scalar 
octet of resonances around 1.4 GeV) are small. They are however somewhat more significant for the 

o 

Tjs {Mn= 545 MeV from the Gell-Mann-Okubo relation). We then favour the standard CHPT scenario, 
where the linear quark mass term is assumed to dominate the pseudoscalar masses, versus the generalized 
one [68]. The M^g that we obtain, around 635 MeV, is larger than the one of the r] physical meson, 
547.45 MeV, that is very close to the bare r] mass calculated from Gell-Mann-Okubo. Our value for M^g 
is very similar to that of ref.[36]! -^rjg = 639 MeV. 

Next, we also show our results for the solution of eq. (l2.34p using the parameters in the S-waves with 
their values determined by fitting experimental data employing always physical masses in the interaction 
kernel. That is, we do not refit them in terms of the bare masses, as done previously. The resulting 
numbers 

M^= 125 ± 4 , Mk= 477 ± 10 , M^g = 633 ± 15 MeV, (3.42) 

are very similar to those already determined. The process of refitting the free parameters in the S- 
waves by distinguishing between physical and bare masses, where it corresponds, just introduces minor 
corrections (compared with the errors quoted) and will not be further considered in what follows. 

Taking into account the expression of the bare masses in terms of the quark ones, eq. (|2.32p . and the 
values in eq. (l3.40p . 



o 2 

rh ° 2 



1 = 27.1 ±2.5. (3.43) 



This number, within errors, is in tli6 bulk of otlicr determinations, v^a — 

25.7 ± 2.6 from 0{p'^) CHPT 

[3j, or with the more refined one of ref. |69j . = 24.4 it 1.5, and also with lattice determinations |30j . 
27.4 lb 0.5. It is worth stressing that our results in this section are predictions without any new free 
parameter. They are a consequence of the strong T-matrices used. 



4 Including 0{p^) CHPT self-energies 

In this section we consider eq. ()2.35p . which reproduces the CHPT self-energies at 0{p^) and the higher 
order contributions are estimated by adding Sp. The r] — rj' mixing is incorporated at 0{p^) by the 
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xlO^ 


m 


[35] 




m 


LI 





~ 0.2 


0.2 ±0.3 [H] 


-0.3 ±0.5 


LI 


0.93 ±0.11 






1.4 ±0.5 


LI 





<o 


0.4 ±0.2 [72] 


-0.2 ±0.3 


Lr 


-0.31 ±0.14 






-0.4 ±0.2 


LI 


0.60 ±0.18 






0.9 ±0.3 


2LI - LI 





~ -0.2 


0.6 ±0.5 


-0.1 ±0.8 


2LI - LI 


0.23 ±0.38 






0.4 ±0.8 


3L7 + LI 


-0.33 ± 0.46 






-0.3 ±0.7 



Table 1: Values for some O(p^) xPT counterterms L[(Mp) x 10^ at the renormalization scale of the p 
mass from the literature. 



counterterm L7 and is fixed to its physical value [3]. Eq. (j2.35p incorporates the combinations of the 
low energy constants ^(5,8) =2-^8 — ^5,-^(4,6) = 2Lg — -L4 and -L(7 g) = SL-r + LZ. eqs. (l2.30p and (12.3ip . In 
table [1] we show values for these 0{p^) LECs at the mass of p resonance, p = 770 MeV. We show in the 
second column the so-called main fit from the 0{p^) fits of ref.[7] and in the fifth one those values from 
the 0{p'^) fits of ref.[70]. In the fit shown of ref.[7] the large Nc suppression of the chiral counterterms L4 
and Lg is used to fix their values to 0, although the precise scale at which this occurs it is not known. The 
error bands in the values of -L4 and Lg in ref.^TOj, fifth column, are intended to cover this uncertainty. In 
this reference the null value for these counterterms is taken at the renormalization scale of the rj meson 
mass. Other fits in ref.[7] are also given where this condition is relaxed and the differences are well 
taken into account by the estimated errors. In the third column we also show the values for L^ and Lg 
from the simultaneous study in ref.^35j of tttt, Kir scattering together with the scalar form factors, the 
masses and decay constants of the lightest pseudoscalars, and K£4^ decays at 0{p^) in CHPT with three 
flavours. The resulting values for such counterterms are found to be small and compatible with zero, in 
agreement with the large Nc suppression. Note that these new values of ref . [35) . with vanishing L4 and 
Lg counterterms, reinforce the results of the main fit of ref. [7], where this was assumed. In refs.|7T| 172] 
chiral sum rules are used to fix the latter counterterms to larger values. These determinations, however, 
suffer of the generally poor (if any) convergence of the three flavour CHPT series [H [3 135] , with sizable 
0{p^) contributions. Because CHPT with three flavours is used only at 0{p'^) the determination of those 
counterterms incorporates large 0{p^) contributions, which are buried in the numbers given. We shall 
comment more on this issue with a precise example below. It should be understood for the third column 
in table [1] that when the value for a given L[ is not shown in the table this is the same as the one in the 
second column. 

We first consider the matching at 0{p'^) of the self-energies given in eqs. ()2.34p and (|2.35p . The vr and 
K self-energies given by eq. (l2.34p at 0(j>^) are 



9 o 2 

Ml = M^{ 



9 ° 2 

Mf^ = Mk { 



32vr2/2 ^^4/i2e 3 327r2/2 ^ 4^2g 

o 2 Q 2 
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with loge = 1. Comparing this expression with the one in eq. (j2.30p one easily concludes that 



where p, is the renormalization scale in CHPT and iJ. was already introduced in eq. (|2.16p . Notice that 
the expressions in eq. (j4.45p are 0{N^), subleading in large N^. This is certainly expected because 
this contribution originates from the loops resummed in eq. (|2.34p . Note that and Lg are 0{N^), 
subleading in the large Nc counting, while L5 and Lg, separately, are of leading order, 0{Nc). However, 
the combination 2Lg — is numerically suppressed, as it is clear from the second column in table [T] from 
ref.^. In ref.jHj one gets the constraint 2Lg — Lg = by requiring that the strange scalar form factors 
vanish in the large Nc limit for s — 00. 
From eq. (j4.45p it follows that 

^(5,8) + 6L(4,6) = . (4.46) 

This is a robust output from our main assumption, namely, the dominance of S-wave rescattering in 
the lightest pseudoscalar self-energies in analogy with the one already well tested for the scalar form 
factors. Note that this relation is independent of the renormalization scheme chosen in the calculation 
of our integrals, contrary to the independent values of 2Lg — Lg and 2Lg — L4 in eq. (|4.45p . E.g. in these 
equations one has the appearance of the number e in the denominators of the logarithms as a remnant 
of our choice. This is removed in eq. (|4.46p . which is also independent of the renormalization scale Ji in 
CHPT. 

One can also extract interesting consequences from eq. (|4.46p . Because of the factor six in front of 
L(-4 g) and since L(5 g) is small in modulus, 

L(4,g) ^ . (4.47) 

This also implies that 

LI ^ hi , (4.48) 

in good agreement with their most recent estimates in ref.[35] from 0{p^) CHPT, third column in table 
[TJ As commented above, the values for L^ and Lg from refs.|7H I72j can suffer from large uncertainties 
due to the large 0{p^) contributions, because only CHPT at 0{p^) was used. As a example, if we take 
the central values for the bare masses in eq. (l3.40p . and evaluate with eq. (l2.30p L(4g) and L(5 g), 

L(4,6) = 0.23 • 10-3 , L(5,g) = . (4.49) 

These numbers clearly violate the constraint in eq. (lOel) . One would obtain L(5_g) ~ -1.3 • IQ-^ if the 
previous value for L(4 g) were used in eq. (l4.46p . This is a clear example that 0{p^) contributions, and 
possibly also higher order ones, are so important in the CHPT series that unless they are taken proper 
and explicitly into account, an estimation of the 0{p'^) counterterms from phenomenology with only 
0{p^) CHPT is severely biased by large higher order contributions buried in the values obtained. 

Next, we apply eq. (l2.35p to study the physical masses of vr, K and rj. As a function of L(4 g) and 

L(5 g), we solve first for the and Mk bare masses in eq. (j2.35p . and fix by the Gell-Mann-Okubo 
relation, eq. (l2.33p . L(7g) is solved from the equation, 

m2 =mI +S^^(Mq) + S^(Mq; m2 ) , (4.50) 
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Figure 6: Left panel: Contour-plot for r„i = nis/rh as a function of i(4.6) and i(5.8)- The square corresponds to 
the second column of table [l] the diamond to the fourth one and the circle to the fifth column. The straight line 
corresponds to the relation in eq. ()4.46p . No values for L5 nor Lg are provided in refs.[7Tl [Tj, and the absence of 
the vertical errorbar in the diamond refers to the fact that ^(5,8) cannot be determined then from these references. 
Right panel: Contour-plot for ^(7,3) as a function of ^(4,6) and i(5,8)- The meaning of the points is the same as 
in the left panel. Our results correspond to the values lying along the straight line within the stripped region, 
indicated by the thick solid line. For more details see the text. 



We use the method of eq. (|3.39p for solving eq. (j2.35p . though the last term on the r.h.s. of eq. (j3.39p is 

o o 

removed since we are considering now the solution for M-k and Mk- In figEt^ we show by the contour 



lines the calculated values for = ms/rh as a function of -^'(4,6) and -^^(5,8) units of 10~^Ej The 
up-left corner where no contour lines are plotted determines a region where no solutions are found. The 
range of values for ^^(4,6) and -^^(5,3) has been chosen to span generously the values of table [T] within errors. 
If one takes into account that ms/rh = 24.4 it 1.5 [69], the preferred region for ^^(4,6) and -^^(5,8) would 
be between the 23 and 26 contour lines. The shadowed areas below and above these lines, moving close 
to the 22 and 27 contour lines, represent the uncertainties of our calculations along them due to the 
variation in the renormalization scale, ^ ~ [0.5, 1.2] GeV, and in the input parameters for the S-waves. 
In the figure, the square represents the values of -^^(4,6) and -^^(5,3) from ref.[7], second column in tabled) 
The full circle corresponds to the values in the last column of the same table and the diamond the value 
of i(4,6) ill the fourth column, refs.[7Tl [72]. Notice that no values for L5 and Lg are provided in these 
references and hence no value for ^(5,3) can be either determined. This is reflected in the absence of the 
vertical errorbar for the diamond. In the figure it is also shown the line corresponding to the constraint 
of eq. (|4.46p . We show in figEja a contour plot for the values of -^^(7,8) that are obtained from eq. (j4.50p . 
If one considers the values L(7^8) = -0-33 ± 0.46 [7j at 0{p^), and L(7^8) = -0-3 ± 0.7 [70j at ©(p^), the 
preferred ^(5,3)) -^(4,6) region would be in the interval ~ [—1,0.4]. Again the shadowed area along these 
lines represents the uncertainties of our calculation. The same points as in figlHK are shown here and also 
the relation of eq. ()4.46p is depicted by the straight line. The stripped area in the two figures represents 
the overlap between the favoured regions on them. Our values lie along the thick solid line defined by 
the intersection of the relation in the eq. ()4.46p with the stripped region. It is certainly remarkable the 
large overlapping between this line and the square point from ref.[7]. Note that our approach and that of 
ref.[7] are completely independent. While refs. [35^17] employ resonance saturation hypothesis to estimate 
the 0{p^) counterterms, which has an important impact on their NNLO results, we have made used of 



*^In the whole manuscript, the values of the LECs displayed in the figures are given in units of 10 ^ and referred to a 
renormalization scale equal to the p mass. 
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the dominant role of rescattering in meson-meson / = 0, 1/2 and 1 scalar dynamics to calculate Sp. In 
addition, while refs.[35l[7] consider simultaneously the lightest pseudoscalar masses and decay constants, 
Ki4 decays, scalar form factors and low energy vrvr, Kit scattering, we have considered all the scattering 
channels in the figlUand the it, K and r/ masses. Taking the intersection of the thick line in figlSh and 
the square from ref.[7j at the one a level, we end with the preferred interval of values 

- 0.15 <L(5,8)-103< 0.35 , (4.51) 

for fi = Mp. Let us recall also our previous estimate for i^(4^6) — in eq. (|4.47p . 

We end this section with a remark. One can also use eq. (l4.45p to derive values independently for 
L(5 8) and i(4,6)- Nevertheless, these values are affected by the choice of the renormalization scheme, so 
that they can differ from those in the standard modified MS scheme of CHPT. In addition, allowing /x 
to change within the range 0.5 — 1.2 GeV the range of values of -^^(5,8) is rather large and does not add 
any real new restriction to what it is already shown in figlH Notice that when we consider eq. ()2.35p the 
0{p^) contribution is given by S^, so that the L[ are considered properly in the same renormalization 
scheme as in CHPT. The same can be said for the relation in eq. (j4.46p as it is invariant under any new 
constants that could be added when removing the infinite part to calculate the integral in eq. (j2.17p . 

5 Comparison with lattice QCD results 




2LI - LI 

Figure 7: Contour-plot for r™ = nis/rh as a function of i(4,6) and i(5,8) obtained by employing only the 0{p'^) 
CHPT self-energies. The point corresponds to the values of ^(5,3) a-nd ^(4.6) from the MILC collaboration [3T| . 

It is worth comparing the result in figE^ with the calculation obtained by considering the self-energies 
only up to 0{p'^) in CHPT, shown in figl?! The difference in ms/rh for the same ^(4,6) and -^^(5,8) is 
typically around 1 — 3 units between both calculations, with larger values when only the O(p^) CHPT 
self-energies are taken into account. A similar trend was also observed in ref.[8] from pure CHPT to 
0{p^), with smaller values as higher orders in the chiral series are included. To O(p^) the values 24.9 
and 24.4 for become, respectively, 24.1 and 23.3 to 0{p^) 
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Two loop diagrams contributing to the self-energies were found very relevant numerically in ref.^, 
larger by several factors than the C(p^) contributions, similarly to what we find. Since diagrams with 
more than one loop were neglected in refs. |31[ [30 ]. worries arise about systematic uncertainties hidden 
in the values determined in these references for low energy parameters. To illustrate this point, we 
show in figlT] by the solid circle the values L(5_8) = (0.16 ± 0.20) • lO'^ and L(4_6) = (0-4 ± 0.4) • 10"^ 
given by the QCD lattice calculation of the MILC Collaboration [31j. The latter reference also predicts 
ms/rh = 27 A ± 0.4. It is worth stressing that this value is very close to the one obtained by considering 
only 0{p^) CHPT, as shown in the figure. Note that 0{p^) CHPT also includes at most diagrams with 
one chiral loop. As a result, we find quite reasonable to expect a reduction in the value of m^/m of one 
unit at least, if higher loop diagrams were included, as we have done. Certainly, further arguments to 
those given in ref.[3T] are necessary to rule out this expectation. 

The MILC Collaboration also provides results for the pseudoscalar meson masses. We present next a 
good fit to these lattice data on the pseudoscalar masses (vr and K), at the level of 1%, making use of 
eq. (j2.35p for the self-energies. We are aware that our parameterization does not take into account the 
taste symmetry violating effects, due to the finite lattice spacing. Nevertheless, we think that our fit is 
valuable since the differences between our full results and those from a fit to lattice data with only O(p^) 
CHPT, eq. (j2.30p . is an estimate of these diagrams with a higher number of loops. It provides then a 
rough quantification of systematic uncertainties that could affect the values for ^(4,6); -^(5,8) ms/rh 
from refs. |30i [31} I32j . It is true that these references include Oijf') contributions at tree level through 
the appropriate chiral counterterms in Staggered CHPT. Nevertheless, since two-loop diagrams, as well 
as one loop diagrams with higher order vertices, are neglected, the procedure is not systematic and some 
quantification is of interest. 

We employ our eq. (j2.35p to reproduce the quark mass dependence of the pseudoscalar masses provided 
by the MILC Collaboration [3T] with three dynamical quarks and equal valence and sea quark masses. The 
systematic procedure to extrapolate these lattice results to the lightest quark masses in the continuum 
is to employing Staggered CHPT [331 [34] with three flavours. Nevertheless, present applications suffer 
of a huge proliferation of free parameters and include only diagrams involving at most one loop, though 
supplied with tree diagrams up to NNLO [31j. 

The strategy is the following. The MILC Collaboration [31j provides pseudoscalar meson masses as a 
function of the bare quark masses used. The current quark masses and the bare ones are proportional 
[30], 

M^=2BoC{a)am , M k= BoC{a)a{m + rus) , M^= -BQC{a)a{rh + 2ms) . (5.52) 

The dependence of the constant C{a) with the lattice spacing can be found in ref.[32] at the two loop 
level, and in ref.[30] calculated up to the one loop order. Previous calculations with C{a) at the one loop 
level |30) gave a value for m*^'^(2GeV) = 76 it 8 MeV, to be compared with the latest value given in 
ref.[32] ^(2GeV) = 87 ± 6 MeV. This substantial shift upwards comes from the improvement in the 
determination of C{a) by moving from the one to the two loop calculation. We use in our fits the two 
loop result for C{a) from ref . [32] . Ref.[31] gives the so called "coarse" lattice runs with dcoarse 

~ 0.12 fm, 

and the "fine" lattices with afim — 0.09 fm, with an error in a about 1.2%. Making use of ref.[32] we 
then determine that C{afine)/C{acoarse) = 2.76/1.85 = 1.49 for n = 2 GeV. From this relation we fix 
C {a fine) in terms of C{acoarse) and take the latter as a free parameter, denoted in the following just as C. 
As the constant C does not depend on the value of the quark mass in very good numerical accuracy (at 
the level of 0.1% [301132] ). we then employ the same value for all the quark masses used in the lattice runs 
with the same a. We solve eq. ()2.35p for the kaon and pion masses given by lattice QCD in terms of the 
bare quark masses used. From the mass equations at the physical point one can obtain i(4,6) and i^{5,8) 
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as a function of and Mk- Hence, the only free parameters in the fit to lattice data are BqC, M-,^ and 

o 2 

Mx- the fit we do not consider the two points of the fine lattice run with nis/m = 0.031/0.0124 = 2.5, 
because it implies a rather heavy pion with a mass of 470 MeV [31] . We then restrict ourselves to the 
case of lighter pions with nis/rh > 5. In figsE^, b, we show the reproduction of the lattice points [3T] and 
of the physical values for the pion and kaon masses, respectively. The size of the squares corresponds to a 
relative error of 1.2%, the one given to a [31]. The circles correspond to our points for the coarse lattices, 
the triangles for the fine ones and the arrows indicate the values of the physical kaon and pion masses. 
It is also worth stressing that the scattering data in figU are simultaneously reproduced!^ From this fit 
we get the values, 

2Lg - = -0.52 ± 0.43 , 2L^ - = -0.20 ± 0.17 . (5.53) 

In figEjwe show by the solid ellipse the one sigma region for -i^(4,6) and -Z>(5,8) obtained in our fit, from 
where one can infer the correlation between these two parameters. We also show the straight line from 
eq. (|4.46p , the square point corresponding to ref. [7] , second column in table [U and by the circle we show 
the values 

2L^ -L^ = 0.16 ±0.20 , 2L^ - = 0.4 ± 0.4 , (5.54) 



obtained by the MILC Collaboration [30^ 131]. The fact that the values in eq. ()5.53p do not lie within the 
favoured region in figl9] is an indication of non- negligible finite lattice spacing effects, already stressed in 
ref.[31j. The MILC values in eq. (j5.54[) are barely compatible, at the level of one sigma, with the thick 
solid line shown inside the stripped region and with the square from ref.[7]- 
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Figure 8: Lattice data from the MILC Collaboration [ST and their reproduction by our parameterization. The 
masses are given in MeV. The size of the squares corresponds to a relative error of 1.2%, the one given to a in 
ref. [31] • The circles refer to our calculations for the coarse lattices and the triangles for the fine ones. The physical 
mass values are signalled out by an arrow. 

We have also performed a fit employing only the self-energies calculated at 0{j,'^) in CHPT. The 
result obtained from the fit to the lattice data is 1L\ - LI = -0.20 • lO'^ and 2Lg - = 0.13 • lO'^ 

In the fit to the unphysical points of lattice the dependence on the bare masses of the pseudoscalar decay constants 
/it, fii and /,,, generally called fp, could play a role. We have then also performed fits to the lattice data recalculating the 
fp in terms of the CHPT expressions for fp as a function of the bare masses, with values for Z/4 and L5 such that 

at the physical point the constants fp have their physical values. As we have checked, these considerations affect little the 
resulting fit and the results are well within the uncertainty bounds quoted before. 
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and these data is quite well reproduced. The difference with respect to eq. ()5.53p is due to the extra 
unitarity loops and higher order effects from the exchanges of tree level resonances. The former effects 
are not systematically considered in ref . [31] . The size of this difference is then an estimate of systematic 
uncertainties that could affect those results from ref.[3lj and that, when added in quadrature with the 
errors given, would enlarge by a factor ~ 2 the errorbar for i^(5,8) eq. (j5.54p . 

rUs 




2LI - L\ 



Figure 9: Contour-plot for = rag/ in as a function of L(4.6) and -^(5 g), showing also the ellipse of the fitted 
values in eq. (|5.53p . The circle corresponds to the lattice extrapolation from Staggered CHPT (31]. The 
rest is the same as in figl6^. 

From the favoured region in figlS^, given by the intersection of the straight line with the stripped 
region and the square point within errorbars, it results 

^ = 22 - 25. (5.55) 

m 

This value is in good agreement with 24.4 it 1.5 from ref.[63. The previous result is somewhat lower 
than the value 27.4 it 0.5 from the MILC Collaboration [30j and ref. [32] . As discussed above, we consider 
that this difference is mostly due to the non inclusion of higher loop diagrams in the Staggered CHPT 
extrapolations used in these references. It follows from this discussion that those calculations are necessary 
to settle this issue. Within the same favoured region of values for ^(4,6) ^i^d -Z>(6,8) ^Iso obtain, 

L(7,8) = -(0.7 -1.1) -10-3 , (5.56) 

and 

M^=116-119, Mii-= 400 -420 MeV. (5.57) 

Finally, we give in table [2] the relative sizes of Y^p and Sp for the same interval of values for i>(5,8)) 
eq. (l4.5ip . and -^^(4,6) = 0. The fact that Sp is larger than Sp- is in agreement with ref.[7j. Indeed, the 
sizes of our self-energies from Sp are rather similar to those determined in this reference to 0{p^). Ref. [8] 
obtains: T,l^/M^ = 0.132 - 0.355, S^/M|. = 0.194 - 0.423 and ^fj^/M^ = 0.234 - 0.521. These values 
are well inside our bulk of results in table [2j Thus, the calculations of ref. [7] to 0{p^), although showing 
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vr 


K 






[-0.010,0.015] 
[0.233, 0.235] 
[-0.05,0.06] 


[-0.01,-0.07] 
[0.352,0.358] 
[-0.03, -0.20] 


[-0.16, -0.18] 
[0.445, 0.460] 
[-0.37, -0.40] 



Table 2: Relative sizes of the 0{p^), 0{p^) and higher order contributions to the self-energies. The 
intervals corresponds to the interval of values for -^^(5,8) eq. (j4.5ip and -Z>(4^6) — 0- 



that this order is much larger than the 0{p^), does not imply necessarily the lack of convergence of the 
chiral series. Our calculation, estimating higher orders corrections by incorporating physical S-waves 
which include both resonant and non-resonant dynamics, gives us values of similar size to those of this 
reference up to 0{p^). In addition, we have to recall the agreement between the straight line in figEh, 
determined from an algebraic matching at O(p^) of eqs. (j2.34p and (j2.35p . and the square from ref.[7j, 
determined numerically from a phenomenological analysis to 0{p^). 

6 Conclusions 

We have undertaken a non-perturbative chiral study of the self-energies of the lightest pseudoscalar 
mesons. In their evaluation the S-wave amplitudes obtained in UCHPT, that reproduce scattering data 
very accurately up to around 1.4 GeV, are employed. These amplitudes generate the exchange of the 
nonet of the lightest scalar resonances, a, /o(980), ao(980) and k, whose contributions to the pseudoscalar 
masses are expected to be relevant since they significantly affect the scalar form factors [26 1 l29 l [73 t [7i]. 
The latter are tightly related to the pseudoscalar self-energies due to the Feynman-Hellman theorem. As 
stressed, their contributions have to be calculated simultaneously with those from the unitarity loops, 
since these resonances due their origin to the large unitarity contributions driven by chiral symmetry. 

We determine that the pseudoscalar masses squared for tt and K are dominated by the bare ones, with 
the self-energy contributions not larger than 10%. For the r] meson they can be larger but still lower than 
around 30%. This favors standard CHPT versus its generalized scenario. We also determine the mass 
of the Tjs, M^g = 635 MeV, in agreement with a previous determination. Next, we take the exact 0{p^) 
CHPT self-energies and resum the higher order contributions within our approach. Matching algebraically 
at 0{p'^) with the calculated self-energies from CHPT one obtains the relation -^^(5,8) — ~6-L(4 g), which is 
independent of the renormalization scale and of our choice of renormalization scheme. Prom this relation 
it follows that ^(4,6) — and thus Lg ~ L\/'^-, in good agreement with the 0{p^) CHPT results from 
ref.[7]. Considering as well the value of the quark mass ratio mg/m = 24.4 it 1.5 from ref.[69], and the 
values for -Z^(7,8) from refs.[71 [TO], we delimit a shorter interval of values for i^(5,8) lying along the line 
-^(5,8) + 6i(4,6) = 0- The value obtained from refs.[35| 17] perfectly overlaps with our determination, and 
the common intersection between these two independents methods gives —0.15 < -^{5,8) " ^ 0.35. This 
implies that we estimate ms/m = 22 — 25 and ^(7,8) = —(0.7 — 1.1) • 10~^. We also warn about higher 
loop diagram effects in the chiral series when evaluating ms/ifi, as they tend to decrease the value for this 
ratio. This is relevant for its recent calculation from lattice QCD [31] as these higher loop effects are not 
systematically accounted for in present extrapolations to physical light quark masses and the continuum. 
Their evaluation is called for. We have also offered a good reproduction of the lattice pseudoscalar masses, 
though the fitted values for -Z>(4,6) and -Z>(5,8) are outside the previous favoured region. This is considered 
as a clear indication of discretization effects not yet negligible in the lattice data, as already stressed in 
ref.[3T]. 
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Other interesting result from our investigation is that although the 0{p^) contribution to the self- 
energies calculated in three flavour CHPT [8j is much larger than the 0{p'^) ones, we still obtain self- 
energies of similar values to those of ref.[8] despite higher order corrections are also included. This fact 
together with our agreement with ref.f7| in the values for -Z>{4,6) and -Z^(5,8); indicate that the SU(3) CHPT 
expansion is not spoiled once the 0{p^) is taken into account. 
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